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Abstract: Certain configurations of D-branes, for example wrong dimensional branes or 
^ • the brane-antibrane system, are unstable to decay. This instability is described by the 

appearance of a tachyonic mode in the spectrum of open strings ending on the brane(s). 
The decay of these unstable systems is described by the rolling of the tachyon field from the 
unstable maximum to the minimum of its potential. We analytically study the dynamics 
■rj" \ of the inhomogeneous tachyon field as it rolls towards the true vacuum of the theory in the 

-^ ' context of several different tachyon effective actions. We find that the vacuum dynamics 

of these theories is remarkably similar and in particular we show that in all cases the 
A ' tachyon field forms caustics where second and higher derivatives of the field blow up. The 

^ ■ formation of caustics signals a pathology in the evolution since each of the effective actions 

considered is not reliable in the vicinity of a caustic. We speculate that the formation of 
caustics is an artifact of truncating the tachyon action, which should contain all orders 
of derivatives acting on the field, to a finite number of derivatives. Finally, we consider 
Cd \ inhomogeneous solutions in p-adic string theory, a toy model of the bosonic tachyon which 

contains derivatives of all orders acting on the field. For a large class of initial conditions 
we conclusively show that the evolution is well behaved in this case. It is unclear if these 
caustics are a genuine prediction of string theory or not. 
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1. Introduction 

Recently effective actions describing the open string tachyon have recieved considerable 
attention in the literature. These effective actions are interesting to study since they 
permit a relatively simple formulation of the complicated open string dynamics in terms 
of classical field theories. In particular, the use of effective actions has triggered significant 
interest in the possible role of the tachyon in cosmology. For example, there have been 
numerous attempts to make use of the tachyon as either the infiaton or as quintessence Q 
(see [@] for reviews of tachyon cosmology). Constraints and shortcomings of these scenarios 
have been discussed in Q. 

Perhaps the most promising application of the tachyon to cosmology comes from D- 
brane inflation |Q]. In the context of D-brane inflation, the annihilation of a brane-antibrane 
pair at the endpoint of inflation is described by the rolling of the tachyon, T, from the 
unstable maximum T = to the degenerate vacua of the theory T = itcxD. The tachyon 
field may roll to the same vacuum everywhere in the space or topological defects may form 
through the Kibble mechanism ||5[. The latter possibility corresponds to the formation of 
lower dimensional branes at points where T = 0. In the regions between the descendant 
branes the tachyon field asymptotically rolls towards the true vacuum of the theory. The 
formation of lower dimensional branes is of some phenomenological importance since the 
resulting cosmic-string-like or higher dimensional defects could be observable remnants of 
infiation Q. In [^, |^] a more radical idea was explored: it was proposed that our own 
observable universe could be such a defect in the higher dimensional spacetime predicted 



by string theory. In this scenario the coupling of the time dependent tachyon condensate 
to gauge particles which will be localized on the descendant brane can provide an efficient 
mechanism for reheating at the end of brane- antibrane inflation. 

Applications of the tachyon to cosmology generically involve the tachyon field rolling 
from the unstable maximum of its potential to the true vacuum of the theory, at least in 
some region of the space, so it is of some interest to study the dynamics in this regime. Often 
only homogeneous tachyon profiles are considered, though a consistent study of the tachyon 
field in cosmology must include spatial inhomogeneities. In the context of the formation 
of lower dimensional branes at the endpoint of D-brane inflation spatial inhomogeneity is 
unavoidable since the field must cross T = at one or more points in the space. However, 
even in the absence of topological defects one expects that spatial inhomogeneities will be 
generated starting from a homogeneous profile by vacuum fluctuations of the tachyon field. 

The full tachyon action contains an infinite number of derivatives acting on the field, 
though the action is not known explicitly to all orders in derivatives. Rather, numerous 
effective actions describing the tachyon have been proposed in the literature on different 
grounds and have been shown to reproduce various nontrivial aspects of the full dynamics 
expected from string theory. Most studies of the tachyon dynamics in terms of effective 
actions have considered actions involving only first derivatives of the field and it is not clear 
how higher derivative corrections may alter the dynamics. Indeed, studies of p-adic string 
theory suggest that there may be fundamental differences between dynamical equations 
with an infinite number of derivatives and those with a large finite number of derivatives. 
In 1^ it was shown that the initial value problem is qualitatively different in these two 
cases. 

In this paper we study analytically the dynamics of the inhomogeneous tachyon field 
as it rolls towards the vacuum in the context of several different effective field theories. 
Our interest is motivated by Q in which inhomogeneous vacuum solutions were discussed 
using the popular tachyon Dirac-Born-Infeld (DBI) action. The authors of [^ found that 
in general the inhomogeneous solutions of the DBI action formed caustics where the second 
and higher derivatives of the tachyon field blow up at some point. In the vicinity of these 
caustics the tachyon DBI action is not reliable since the higher derivative corrections which 
have been ignored will become important. It is thus not known if the formation of caustics 
is a genuine prediction of string theory or an artifact of the DBI action for the tachyon. To 
shed some light on this question it is interesting to study inhomogeneous vacuum solutions 
of different tachyon effective actions and look for caustics or similar singularities. 

The finite time formation of caustics found in Q must be distinguished from the finite 
time divergences in the first derivatives of the tachyon field at points where T = [0, 10 1 



which corresponds to the formation of topological defects. In both cases the divergences in 
the derivatives of the tachyon field lead to divergences in the energy density. In the case 
of topological defect formation, however, this finite time divergence in the energy density 
has the form of a delta function and hence leads to finite total energy [0]. We point out 



that similar behaviour was found in [11| where the decaying tachyon on the supergravity 
background corresponding to space-like branes in string theory was studied. In this case 
the first derivative of the tachyon field and the energy density of the tachyon matter both 



divergence in finite time producing a space-like curvature singularity in the metric. 

The organization of this paper is as follows. In section |2| we consider three tachyon 
effective actions, all of which contain only first derivatives of the field and show that all 
three display surprisingly similar dynamics near the true vacuum of the theory. In section 
m we briefly review how this universal vacuum structure leads to the formation of caustics. 
In section Q we consider the vacuum dynamics of the tachyon in the context of yet another 
effective action which contains terms with second derivatives of the field and show that the 
problem of caustic formation is not ameliorated. In section |5| we briefly review the bosonic 
tachyon action of p-adic string theory and study small spatial inhomogeneities about a time 
dependent solution which rolls towards the minimum of the potential (which is unbounded 
from below). We find that for a large class of initial data our perturbative expansion is 
reliable throughout the evolution and the solutions are well behaved. Finally we conclude 
and discuss the difficulty of interpreting caustic formation physically. 

2. First Derivative Tachyon Effective Actions and the Eikonal Equation 

For simplicity we restrict ourselves to studying real tachyon fields in 1 + 1-dimensional 
Minkowski space with metric 'q^^dx^dx'^ = —dt^ + dx'^. In the context of brane annihilation 
the rolling of the tachyon in 1 + 1 dimensions from the unstable maximum to the ground 
state describes the decay of an unstable Dp-brane to either the vacuum (in the case that 
the tachyon rolls to the same vacuum everywhere) or to a brane of codimension one (in 
the case that kinks form). We work in units where a' = 2. 

2.1 The Tachyon Dirac-Born-Infeld Action 

We first consider the tachyon Dirac-Born-Infeld action |12, ^, 14 1 



£ = -U{T)^l + d^,Tdt'T (2.1) 

with potential 

where Tp is the Dp-brane tension. Solutions of this theory have been widely studied in the 
literature |^, 0, |9|, |l^, 15, 16]. This action can be obtained from string theory in some limit 



| 17|| and has been shown to reproduce numerous nontrivial aspects of the full string theory 
dynamics [0, IH. 

The action ( ^.l] ) should be thought of as describing tachyon profiles which are "close" 



^ to the homogeneous rolling tachyon solution 

T{t) = 2 sinh^i (t+c*/^ + r^e"*/^) . (2.3) 



We are interested in the dynamics of ( |2.1| ) close to the vacuum T — > itcxD, U{T) -^ 0. 
The vacuum structure of this theory has been studied analytically in lO, O] which we 
briefly review. The vacuum structure is most easily described in the Hamiltonian formalism 



""^That is to say taking T+, T- to be slowly varying functions of x'^ . 



since the Lagrangian does not survive the hmit U{T) — > 0, but the Hamiltonian does. 
Defining the momentum conjugate to T as 11 = dC/dT ^ the Hamiltonian is given by 

H = f (fxTi, n = yjip + c/(r) Vi + r'2. 

In the hmit as U{T) — > one finds 

n = nVi + r'2. (2.4) 

In the vacuum Hamilton's equations of motion are 



t = Vl + T'2 (2.5) 

and 

n = 5. (n^^) (2.6) 

for H > 0. Equation ( |2.5| ) is precisely the eikonal equation 

a^ra^T + 1 = o. (2.7) 



It was observed numerically in [^ that the general solutions of ( p.l[) tend towards the first 
order equation (^). 

2.2 The Boundary String Field Theory Action 

The action 

C = -V{T)F (a^TS'^T) (2.8) 

where the potential is 

y(r) = V2rpexp('-^') (2.9) 

and the kinetic term is 

F(.) = ^^ (2.10) 

can be derived from boundary string field theory (BSFT) assuming a linear profile T = 
a + Ufj_x^, and summed to all orders in u^ |jl^. This action should therefore be considered 
reliable only when describing profiles where second and higher derivatives of the tachyon 
field are small. 

The dynamics of the action (|2.8D have been discussed in pO] which we briefly review. 



Before attempting to describe the vacuum dynamics of the theory (2.8) we discuss some 
limiting behaviour of the function F{z). 



^Here and throughout this paper T = doT = dtT and T' = diT = d^T. 



At small z the function F(^z) has a Taylor expansion 

F{z) ^ 1 + ln(4)z + • • • 
and at large positive z the function Fi^z") has the behaviour 



F(z^ 



TTZ 



as z 



oo. 



(2.11) 



It is also noteworthy that F{z) is singular at z = —n (n = 1, 2, • • •). Of particular interest to 
the ensuing analysis will be the limiting behaviour of F{z) near the singular point z = —1 



F{z) 



2(^ + 1)' 



F'{z) 



1 



2{z + l) 



F"{z) 



-1 



+ 1)3- 



(2.12) 



We now proceed to study the dynamics of the action (2.8). The equation of motion 
which follows from (12.81) is 



d''TduzF"(z) + d'^duTF'iz) 



^TzF'{z) + ^TFiz) 



d^Td^T. We will assume initial conditions such that z 



= (2.13) 

initially. Further- 



where z 

more, we assume that z does not cross the singularity at z = —1 ^. 

The homogeneous solutions of (2.13|) have the asymptotic behaviour T ^^ =bt at late 



times 1 21] so we expect that the tachyon will roll towards the vacuum T -^ ±00 as t — > 00. 
We consider first the case where z does not approach —1 as T — > 00. In this case F{z), 
F'{z) and F"{z) are well-behaved and the leading contribution to ( 2.13|) is 



F{z) - 2zF'{z) = 0. 

It is easy to show that the quantity F{z) — 2zF'{z) is greater than or equal to for all 
z > —1 and that F{z) — 2zF'{z) ^ as z ^ -|-cxd (see equation ( 2.11 )). We conclude 
then that for solutions where z is increasing the equation of motion requires z — > -|-cxd as 
T ^±00. 

We now consider the possibility that z ^ — 1 at late times. Near z = —1 the function 



F{z) and its derivatives are given by (2.12) and hence the leading contribution to ( p. 13 ) is 

1 



d^Td^zF"{z) - -TzF'iz) = 0. 



-1. Since z -|- 1 —> as T 



(2.14) 
±00 it is 



This equation can be satisfied for profiles where z 

reasonable to assume that we can write z + 1 = f{T). This is consistent with (2.14) which 

implies 

z-hi~y(r)^/2^ 

It is important to stress that, barring any asumptions on the form of the solution, the 
equation of motion ( p. 13 ) requires that d^Td^T tends to either -|-oo or —1 as T ^ itcxD. 
The former case describes tachyon kinks of the form T{x) = x/e to be understood in the 



^This is not a particularly restrictive assumption since if z did cross the singularity it would lead to 
infinite action. Note that this does not exclude the possibility that z —> —1 as V^(r) —> 0. 



limit that e — > 0. The latter case describes the rolling tachyon and corresponds to the 
vacuum structure we are interested in. We conclude then that close to the vacuum the 
tachyon field obeys the eikonal equation 

df^Td^T + 1 = 0. 

2.3 The Lambert and Sachs Action 

The action 

C = -V{T)K{d^Tdf'T) (2.15) 

where the potential is given by (^]^) and the kinetic term is 

K{z)=e-^ + V^eif{V^) (2.16) 

was proposed in p^ . This action is derived by requiring that the profile 

I [X) = xsm 



be an exact solution to the equations of motion. The residual freedom is fixed by demanding 
that £ = v2Tpexp(— T^/4) for constant profiles for agreement with (|2.8| ). 

Before proceeding to study the dynamics of (2.15) we consider the asymptotic be- 



haviour of the kinetic term ( ^.16| ). For z > the function K{z) is remarkably similar 
to the BSFT function F{z) given by ( p. 10 ). Near z = the function K{z) has a Taylor 



expansion 

K{z)^l + z + ---. 
For large positive values of z the function K(z) has the aymptotic behaviour 

K{z) -. V^z, K\z) -. ly^, K"{z) -. -^ (2.17) 

so that K{z) tends to infinity slowly as z — > +oo while K'(z), K"{z) tend to zero. Finally 
we consider the behaviour of K{z) for large negative z. Using the formula 

^erf(iy) ^ ie^' (^ + ^ + ^(y"')) 
for y — > cxD it is easy to show that 

p\z\ p\z\ \z\ 

K(z) -^ ----, K'iz) -^ +^, K"{z) -^ -^ (2.18) 

^ ' 2|z|' ^ ' 2|z|' ^ ' 2\z\ ^ ' 

for large negative z. Note that though K{z) has very similar behaviour to F{z) and ^/\ + z 
for large positive z, these three kinetic terms have very different behaviour on z < 0. 



The equation of motion which follows from ( |2.15| ) is 

hrzK'{z) + \ 



df'Td^z K"{z) + df'd^T K'{z) - -TzK'iz) + -TK{z) = (2.19) 



where z = d^Td^T. By construction equation (2.19) has the homogeneous solution 



T{t) = r+e*/2 + r_e"*/2_ 

Comparing this to the homogeneous solution (^]^) demonstrates that a field redefinition is 
necessary to make accurate comparisons between the two theories. Notice that in the case 
of the homogeneous solution the tachyon field rolls towards the vacuum T — > ±00 at late 
times with z — > — 00. 

We proceed now to study the vacuum structure of the theory ( 2.15|) . First we consider 



the possibility that z tends to some finite positive value as T ^ ±00. In this case the 
leading contribution to ( ^.19 ) is 



TzK'{z) = \tK[z) (2.20) 



or 



This suggests that we require z — > +00 if 2; > and indeed, using the limiting behaviour 
( 2.17 ) the leading contribution to ( 2.19| ) for T ^ itcx) with z -^ +00 is satisfied. 



We consider now possibility that z approaches some finite negative value as T — > ±00. 
Assuming z 7^ — cxd then K{z), K'{z) and K"{z) are all well behaved and the leading 
contribution to ( 2.19| ) is still ( 2.20| ) or equivalently = e'^' which has no solution. 



Finally we consider the possibility that 2 — > — cx) as T ^ ±00. In this case the limiting 
behaviour ( 2.1§| ) is applicable and K{z) ~ —K'{z) ~ K"{z) so that the leading contribution 
to ( p^ ) is 

2d^Td^'z + Tz = 0. (2.21) 

Since z ^ — 00 as T — > ±00 it is reasonable to assume z = f{T) in this regime. This 



assumption is consistent with ( 2.21 ) and we find /(T) = — T^ up to an arbitrary additive 
constant which is irrelevant as T^ ^ 00. One may verify that given these assumptions 
the leading behaviour ( 2.21| ) is order T^ while the terms we have disregarded in ( 2.1S| ) are 



order T so that this series of approximations is self-consistent. We conclude that for z < 
the equation of motion requires that the tachyon field obeys 

^f,T^^'T + T^ = 0. 

Defining T = exp(r) we find that T obeys the eikonal equation 

d^fdf'f + 1 = 0. 

We stress that the equation of motion ( |2.19|) requires either z -^ +00 or 2; — > — cx) as 
T — > ±00. The former case describes tachyon kink solutions while the latter case describes 
the rolling tachyon and corresponds to the vacuum structure we are interested in. 



3. Caustic Formation 



We have found that ah three actions (^,2.S,2.15|) predict that the tachyon field is described 



by the eikonal equation (pj) close to the vacuum. Inhomogeneous solutions of ( |2.7| ) were 
found in P] for arbitrary Cauchy data by the method of characteristics. The solution is 
given along a set of characteristic curves which, in 1 + 1 dimensions, are defined by 

x{q,t) = q - ^M^^t (3.1) 

^^' ' ^l + Tl{qY ^ ' 

where the parameter q defines the initial position of the curve on the x axis such that 
x{q,t = 0) = q and Ti{x) is the Cauchy data at time t = 0. The parameter q should be 
thought of as labelling the curves. The value of the tachyon field along a given characteristic 
curve is 

It is also worth noting that the derivatives of the field are constant along the curves 
T'{q,t) = Tl{q) and t{q,t) = ^l + Tl{qY. From (U) it is clear that in general for initial 
data where T/' ^ there will be curves originating from different initial points on the x 
axis which will cross in some finite time. At points where the characteristic curves intersect 
the field becomes multi-valued since evolving ( |3.2| ) along two different curves leads to two 
different values of T at the point of intersection. This corresponds to the formation of 
caustics and signals a pathology in the evolution. To illustrate the problem, consider T" 
along a characteristic curve 

T"{q,t) ^"^^^ 



1 n'fa) / 



At some time tc and point qc where Tl'(qc) > the denominator vanishes and the second 
derivative blows up. For a given qc the caustics form at time 

, il + Tl{qc)y/' 
" ~ r/'(ge) • 



In each of the three effective actions considered above (p?l,2.8,^.15) this caustic formation 



signals a breakdown of the theory since terms in the action involving second and higher 
derivatives of the field which have been neglected will become important. 

4. Vacuum Dynamics in a Higher Derivative Action 

In light of the analysis of the preceeding section it is tempting to speculate that the caus- 
tic formation described above is an artifact of the derivative truncation which leads to 
(2.1, 2.8 , 2.15| ) since these effective descriptions break down in the vicinity of a caustic. 



We would like to reconsider the effective description of the tachyon near the vacuum in 
the context of an action which does contain terms with second derivatives of the field. 



One might consider attempting to generalize the action (|2.8| ) to a profile of the form 
T = a + u^x^ + Vfj_yX^x^ to obtain an action which is valid for profiles where d^T ^ and 
dfidiyT 7^ without constraint on the size of the first and second derivatives. However, 
deriving such an action (following [jl9[) would involve performing a path integral which is 
no longer Gaussian. On the other hand, the superstring tachyon effective action has been 
calculated in an expansion in small momenta (derivatives) around a constant profile up to 
six orders in derivatives {0{d^)) in |23]. The action, in units where a' = 2 and truncated 
to fourth order in derivatives, is 

C = -V{T) ( 1 + ln(4)a^ra^T + (i(ln(4))4 - C(2)) {d^Td'^Ty + C{2)Td^duTdi'Td''T 

+ [(2C(2) - 8(ln(2))2) + ^r2J d^d^Td^'d'^T ) (4.1) 

where ("(z) is the Riemann zeta function. The action (4.1) is valid to all orders in T but 



for small derivatives and hence naively does not seem applicable to studying the vacuum 
dynamics where, the analysis of the preceeding sections would suggest, d^Td^T is order 
unity or larger. However, one might expect the vacuum structure of the theory to be 
significantly modified due to the presence of terms like T^d^dyTd^d'^T in the action ( |4.l] ) 
which are vanishing on linear profiles but are large in the vacuum for nonlinear profiles. One 
might hope that the vacuum structure of the theory (^]^) is such that when one considers 
Cauchy data where the derivatives of T are small but nonzero, then the derivatives stay 
small throughout the evolution and the approximations which lead to ( [4.1[ ) remain self 
consistent. This turns out not to be the case, as we will show. 

The equation of motion which follows from ( [4.1D is cumbersome and we do not write 
it out here. In the vacuum, as T — > itcxD, assuming that the derivatives of the field are well 
behaved ^ the leading contribution to the equation of motion is 

d^'d^Tdf^TdyT = (4.2) 

where the subleading terms are of order T^^ and smaller. 
We note that (|4.2| ) may be re-written as 

d^Td^" (d^Td^T) = 

so that the solution set of ( [4. 21) contains the solutions of the first order equation 

duTd^T = K 

where k is a constant which may be set to +1, —1 or by rescaling T. The case k = —1 
corresponds to the eikonal equation and we conclude that ( |4.2D does admit solutions with 
caustic formation in the case of Cauchy data such that Tj(x)'^ = 1 + T/(a;)^. However, 
equation ( [4. 21) is second order and of course it is not necessary to constrain the Cauchy 
data in this manner. For completeness we discuss the construction of more general solutions 
of ( [4. 2]) in the appendix. 



^If we drop this assumption then ([4.l|) ceases to be a reasonable description of the dynamics. 



It is beyond the scope of this paper to perform a comprehensive analysis of the deriva- 
tive singularities admitted by ( |4.2| ). We have shown an explicit reduction of (^^) to the 
eikonal equation, which exhibits caustics. As discussed in the appendix, we believe more 
general solutions of ( [4.2| ) exhibit similar pathologies. For our purposes, this is sufficient 
to demonstrate that the higher derivative action ( |4.1| ) does not ameliorate the problem of 
caustic formation. In light of these results, it is tempting to speculate that it is necessary 
to use an action which contains all orders of derivatives acting on the field to obtain a fully 
consistent description of the dynamics of the tachyon near the vacuum. 

As an illustrative example we will discuss inhomogeneous solutions in a p-adic string 
theory, a toy theory of the bosonic string tachyon which contains an infinite number of 
derivatives acting on the field which is known explicitly to all orders in derivatives. 



5. Inhomogeneous Solutions in p-adic String Theory 

The action of p-adic string theory is |24] 



S 



p^ 



g'^ip- 1) 



d^x 



,p-h9^9^^(j) + 



1 



0+1 



(5.1) 



P+l 

where cp is the open string tachyon, g is the open string coupling constant and, though the 
action was derived for p a prime number, it appears that p can be continued to any positive 



integer (the action makes sense even in the limit p ^ 1 |25[| ). The differential operator 
p-^OfiO ^g ^Q |-|g understood as the series expansion 



P 



-^d^d'' 



e^p(-^ln{p)d,dA=f2 



1 Y 1 



The equation of motion which follows from ( p.lD is 



P 



ji^ML 



and the potential is 



V{<P) 



g^{p- 1) 



p+l 



AP+l 



y {d,dn' 



(5.2) 
(5.3) 




>• 0.5 




Figure 1: Plot of the potential V{4>) for Figure 2: Plot of the potential F(0) for odd 
even p. p. 
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The cases of odd and even p are qualitatively different. For odd p the potential is an 
even function of (/> and for even p the potential is an odd function of (p. In both cases the 
potential is unbounded from below. In the case of even p the perturbative vacuum is at 
(j) = 1 and in the case of odd p there is an equivalent false vacuum at cj) = —1. In both 
cases the true vacuum of the theory is at = 0. Figures [l| and |2| show plots of the potential 
( |5.3| ) for the cases of even and odd p respectively. 

We note that the action ( |5.1| ) is intended to describe the bosonic string tachyon, 
rather than the superstring tachyon which we have considered in our previous analysis 
(|2Jj2.8| , 2.15 .4l|). Furthermore, (BJ.) is a toy model which is only expected to qualitatively 
reproduce some aspects of a more realistic theory. These points should be kept in mind 
when comparing analysis using (|5.1| ) to the previous analysis using ( ^T , p.8 , p.l5| , 4.1 ). That 
being said we note that there are several nontrivial qualitative similarities between p-adic 
string theory and tachyon matter. For example, near the true vacuum of the theory (p = 
the field naively has no dynamics since its mass squared goes to infinity ^. This is the p-adic 
version of the statement that there are no open string excitations at the tachyon vacuum. 
A second similarity between p-adic string theory and tachyon matter is the existence of 
lump-like soliton solutions representing p-adic D-branes [^]. The theory of small fiuctua- 
tions about these lump solutions has a spectrum of equally spaced masses squared for the 



modes j27], as in the case of normal bosonic string theory. On the other hand, there are 
some important differences between the theory {5A) and (pri, plg , 2.15 ,4.1). In the case of 
tachyon matter the vacuum is at infinity and the tachyon never reaches this point, whereas 
in the case of the p-adic string the vacuum is at a finite point in the field configuration space 
and homogeneous solutions rolling towards the vacuum typically pass this point without 
difficulty m. In fact, the numerical studies of |8[ found no homogeneous solutions which 
appeared to correspond to tachyon matter (vanishing pressure at late times). 

Keeping in mind that the connection between (5.1) and ( ^T , p.8 , 2.15| ,4.1) is not entirely 



clear we proceed to study inhomogeneous solutions of (5.1) as an example of a theory with 
an infinite number of derivatives which may have some qualitative similarities to the string 
theory tachyon. At first glance it is not immediately clear how to proceed since ( |5.2D is 
difficult to solve for profiles with nontrivial dependence on more than one variable. The 
simplest solution is to study small fiuctuations about some known time-dependent solution 
of ( [5.21) . We are interested in the dynamics near the vacuum and [^ found anharmonic 
homogenous oscillations near = so one might consider small spatial inhomogeneities 
about these solutions. However, these anharmonic oscillations cannot be found by solving 
the linearized equation of motion and hence we should not expect to be able to study small 
inhomogeneities near the vacuum by linearizing about these oscillators. One might consider 
the closest analogy to the solutions found in the preceeding sections to be small fiuctuations 
about a time dependent solution which interpolates between the unstable vacuum cp = 1 
(or also </> = — 1 in the case of odd p) and the stable vacuum i;^ = 0. However, it was shown 
in [S, 



that no such time dependent solution to (5^) exists . We therefore will consider 



^Reference [p| found anharmonic oscillations around the vacuum by numerically solving the full nonlinear 
equation of motion. However, these solutions do not correspond to conventional physical states. 

^More precisely, it was shown that there exists no homogeneous nonnegative bounded continuous solution 
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inhomogeneous fluctuations about the rapidly increasing solution 



Zpmp 



(5.4) 



We stress that this solution does not roll from the unstable maximum of the potential to 
the true vacuum of the theory (indeed no such solution appears to exist) and thus the 
ensuing analysis may be of limited relevance since the connection between p-adic string 
theory and tachyon matter is unclear. The homogeneous solution (|5.4| ) does bear some 
qualitative similarity to the solutions of (2.1,p.8, 2.15 ,4.1) in the sense that this solution 
represents the tachyon rolling down the potential though in the case of (^^) the tachyon 
rolls towards y — > — oo and not ^ = 0. 



Writing (j)(t, x) = (j)o{t)+5(j){t, x) with (j)o{t) given by (5.4) the equation of motion (| 



IS 



to linear order in 



where 



•Q. The particular solutions of (^) may be written as 

H\{t,x) = (j)o{t)Kx{iat)x\(x) 



(5.5) 
(5.6) 



Q 



'lp2 



1 



2 plnp 
and K\{z) is any solution of the Hermite equation 



dz 



dz 



0. 



(5.7) 



Plugging (|5.6| ) into ( ^.5| ) one finds that the spatial modes xxi^) ^^^ determined by the 
equation 

P~^^^X\ix)=P^'^Xxix)- (5.8) 

Before attempting to solve (|5.8D for the spatial dependence of the particular solutions 



of (|5.5| ) some comments are in order concerning the separation of variables (5^). It is 



straightforward to show that the solutions of ( |5.5D separate as (5^) using the identity |27| 



=«9? 



Kx{iat)e^ 



(1 - 4a6)-^/2 ( 1 ^ 



Aaa^ \ '"^ 



X Ky. 



iat 



^{l-4ab){l-4ab + 4aa'^) 



exp 



bt^ 



l-4ab 



(5.9) 



This identity was considered in |2^ for A = n = 0, 1,2,--- and Kn{z) = Hn{z), the 
Hermite polynomials, though it holds for arbitrary A which is most easily shown by noting 
that the two solutions of (|5.7D, Kx{z) = Px{z) and Kx{z) = Qx{z), have contour integral 
representations |^ 



Ci 



of (5.2) with (j){t -> -oo) = 1 and (j){t -^ +oo) = 0. 
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and 

Qx(z)= [ ^r^-ie-«^+2«^ (5.11) 

Jc2 2vri 

where the curves Ci and C2 are given in figure y. The standard definition of the Hermite 
polynomials for A = n = 0,l,2,--- is given by 

Hn{z) =n\{Pn{z)+Qn{z)) 

,1 f ^f:-n~l-e+2iz 



= n\ / ^:^^-"-ie-* +^*^ (5.12) 

where C is the contour given in figure Q. For subsequent calculations we will be particularly 
interested in defining the Hermite functions H\{z) on A < 1. For Re(A) < we can contract 
the path of integration C to the origin and write 

^^('^ = f(Zxjl dCr'-'e-^'-'^^ (5.13) 



where the normalization has been chosen to agree with |^0[. Note that when Re(A) > —1 
we can also represent the Hermite functions by a real integral [^] 

H,{z) = ^l^ e"^"e cos (2.^ - ^) (5.14) 



which coincides with the standard definition of the Hermite polynomials when A = n 
0, 1, 2, • • •. With these conventions the Hermite functions are normalized so that 

_ 2^r(l/2) rrl(r,._tltlB. 



.^v. 



e 



^\j- 



Figure 3: The contours of integration Ci Figure 4: The contour of integration C . 

and C2. 



We now proceed to determine the spatial dependence of the modes (5^). Equation 
|5.8|) has solutions 

Xx{x) = ax cos(u;Aa;) + bx sin(a;Aa::) (5.15) 



where ux = y^2(1 — A). For initial data which are periodic on some interval [—L, +L] 
one takes ux = vrm/L where m = 0, 1,2, • • •. With this choice the degree of the Hermite 
function is 



ir'^m? 



' = ' 2L^ 
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For the zero mode m = we have A = 1 and for ah other m we have A 7^ 0, 1, 2, • • • unless 
L is chosen so that V2L/7: is an integer. 

The Hermite equation ( |5.7D is second order and we expect to be able to find two linearly 
independent solutions for each A. It is most convenient to choose one of these to be the 
Hermite functions as defined by ( |5.12 ,5.13, 5.14| ). To obtain a second solution we note the 



Wronskian formula [30| 



W[H,{iz),H,i-iz)] = =f^e'^\ (5.16) 

so that 

6<p[+\t,x) = Mt)Hxi+iat)x[^\x), 
6(p[-\t,x) = Mt)Hx(-iat)xt\x) 

are linearly independent for A 7^ 0, 1, 2, • • •. The linear independence of these two solutions 
fails in the exceptional case A = 1 for all L. The linear independence of these two solutions 
may also fail on other values of m if ^/2L/7: is an integer. For simplicity we exclude the 
case of integer \/2L/7: from the present analysis and take A = 1 to be the only special 
case. (We note, however, that is is straighforward to extend our analysis to include integer 
values of \'zL/7r.) In the case A = 1 (m = 0) we can take, as one solution of (5.7), the 
Hermite polynomial ( 5.12| ) 

Hi{z) =2z. 

The second solution may be written formally as 

Piiz)-Qiiz) 

where the functions P\ and Q\ are given by ( ^.10| , ^rTl| ). The function obtained in this 
manner does not have a simple closed form expression as does Hi(z). 

Putting the results of this section together we write the general solutions of (|5.5D as 



(5(/>(t,x) = Mt) Re [Pi{iat) - Qi{iat)] i aj,^^ + 0o(i) | 4^^ (5.17) 

Mt) E^=l H['l,my2L^ (t) X^^ (X) + Mt) E™=1 ^J1.„V2L^ (*) ^^^ (^) 

where we have defined the functions 
and 



^(2),,,_ r(-A/2) 1 

2^r(-l/2) 2m 



^A Hi) = oAw .^ — (Hxiiat) - Hxi-iat)) . 



The normalizations have been chosen so that Hj^ (0) = 1 and dtH^^ (t)|t=o = 1- The 
spatial modes {xm} are 

X^)(x) = a« cos (^) + 6« sin (^) 
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for i = 1,2. The coefficients {am , bm} determine the fourier expansion of Scp/cpo at t = 

(2) (2) 

and the coefficents {am \ bm } determine the fourier expansion of dt 



at i = 0. 



For practical numerical computations it is simplest to consider initial data where Og = 
to eliminate the function Pi{iat) — Qi{iat) from the expression (5.17). For example, this 
will be true for any initial data such that 5<j)/(j)Q at t = is an odd function of x. 

Recall that the linearized equation ( ^ ) is valid for |5(/)(t, x)/(/)o(t)| <^ 1. The stability 
of the solution (|5.17D therefore depends on the asymptotic behaviour of the functions 



r(l). 



r(2) 



H^ (t) and H^ (t) at late times. It is straightforward to show that H^ (t) and H^ (t) 
have the asymptotic behaviour 



r(2)/ 



r(i) 



Hi'>{t)^Hi^>{t)^t 



(2), 



for large t. That these two solutions become linearly dependent as t 
anticipated from (5.16). We find, then, that for modes with A = 1 



k(i) 



cxD was to be 
^ < one 



has 5(j)m{t,x)/(j)o{t) — > and the linearized approximation is stable. On the other hand, 

> are increasing functions of time and hence the linearized 



modes with A 



'217' 



,(1) 



.(2) 







approximation will eventually break down. In the case that L < 7r/v2 and 
then all of the modes are decreasing and the perturbative expansion is reliable throughout 
the evolution. In this case we can make definitive statements about the absence of caustics 
or similar singularities in the theory ( |5.1] ). In fact, for initial data which satisfy these 
restrictions, the profile (f){t, x) becomes homogeneous at late times. One might argue that 
(15.17) implies the absence of caustic formation even for more general initial data (for 
example when L is large and many of the modes 64>m{t,x)/cj)Q{t) are increasing) since in 
this case the homogeneous zero mode of 5(j)/(j)Q increases faster than all other modes. 



0.001 - 



o 



-e- 



■0.001 




0.005 




Figure 5: Plot of the spatial profiles o<pi (po 
for a series of increasing time steps. The ini- 
tial data is such that the linearized approxi- 
mation is valid throughout the evolution. 



Figure 6: Plot of the spatial profiles oqy/cpo 
for a series of increasing time steps. The ini- 
tial data is such that the linearized approxi- 
mation breaks down at late times. 



Figure ^ shows a plot of the spatial profiles d(j)/(j)o given by ( |5.17| ) for a series of 
increasing time steps for the case p = 2, L = 1. The initial data are 6(I){0,x)/(J)q{0) = 
^Recall that we are excluding the case A = 0. 
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COlxe-^^/O-"^)" and dti6(f>{t,x)/(j)oit))\t=o = 0.001 sin(7rx/L) on -L < x < +L. For this 
choice of initial data a^ = Oq ~ ^ ^^^ A < for all m. The perturbative expansion is 
stable for this example and the field tends towards a homogeneous profile at late times. 
Figure ^ shows the same plot for the initial data 6(l){0,x)/(j)o{0) = 0.01 x e^*-^'^'^' and 
dt{5(f){t, x) / (j)Q{t))\t=o = 0.01 on —L < x < +L. With this choice of initial data 



(2) 

at late times since Og 7^ and the linearized approximation will eventually break down. 
Finally we comment on the interpretation of the Cauchy problem for the linearized 



differential equation (5.5). We have found a general solution of (p.5|) for which we are 



free to specify the initial field 6(f)/(j)Q at i = and the time derivative dt{d(p/(l)o) at t = 0. 



This may seem surprising since (5^) contains an infinite number of time derivatives and 
one might naively expect to have the freedom to specify an infinite number of initial data 
d^ {5(/)/(pQ)\t=o for n = 0, 1,2, •••. However, the initial value problem for homogeneous 
solutions of (|5.2| ) was studied in [^ and it was found that the equation of motion itself 
imposes an infinite number of consistency conditions on the initial data which one can 
consider. In fact, it was speculated in [^ that the space of allowable initial conditions for 
( |5.2| ) may be finite. This conjecture seems consistent with our results. It is interesting to 
note that ( |5.5D seems to be an example of an equation containing an infinite number of 
derivatives whose solution space is surprisingly similar to that of equations containing only 
two time derivatives. 

6. Conclusions 

We have studied several different effective actions describing the open superstring tachyon. 
In the case of actions containing only first derivatives of the tachyon field we have studied 
three different theories proposed on different grounds. These three actions are expected to 
be valid in different limits of string theory, however, we find that that vacuum structure of 
all three theories is remarkably similar. In particular, we have shown that all three theories 
lead to the formation of caustics where the field becomes multi- valued and where second and 
higher derivatives of the field blow up. Each of these three actions cannot be trusted in the 
vicinity of a caustic since higher derivative corrections which have been neglected become 
important. We considered also an effective action containing second order derivatives of 
the field and found a similar structure of derivative singularities. Finally, in the context 
of p-adic string theory, we studied small inhomogeneities about a time dependent solution 
which rolls from the unstable vacuum to infinity in field configuration space. In this case 
we found that for a broad class of initial data the linearized approximation is reliable 
throughout the evolution and we could conclusively show the absence of caustics or similar 
pathologies. This result seems suggestive that the formation of caustics is an artifact of 
truncating a theory with an infinite number of derivatives, however, the connection between 
p-adic string theory and tachyon matter is still unclear. 

Since we have restricted our analysis to the study of effective actions we cannot con- 
clude if the phenomenon of caustic formation is a genuine prediction of string theory or 
not. If this prediction is borne out by string theory it will be necessary to find some 
physical interpretation of the caustics, though none is obvious to us. As pointed out in B 
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such an interpretation could depend on the dimensionahty of the theory. If these caustics 
are in fact a real prediction of string theory it will also be necessary to find some way to 
predict the field value in the vicinity of a caustic. The situation is qualitatively similar 
to the formation of shock waves in nonlinear gas dynamics. In that case the multi-valued 
continuous solution in the vicinity of a shock is replaced with a single valued discontinuous 
solution using the Rankine-Hugoniot jump condition to quantify the discontinuity in the 
field. It is possible that some similar auxiliary condition could be used to predict the value 
of the tachyon field in the vicinity of a caustic. 
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APPENDIX: The Legendre Transformation 

In 1 + 1 dimensions equation ( |4.2| ) takes the form 

f2f ^ j,i2j,ii _ 2fT'f' = (A-1) 



To construct general solutions of ( |A-1| ) we preform the Legendre transformation (see |31| 
for example) defined by the relations 

e = r', 77 = t, T{t,x)+uj{ii,i)=xi + tn (A-2) 

where u){r], ^) is a new dependent variable and {77, ^} are the new independent variables. It 

follows from ([A-2| ) that 

duo duj , , , 

It is straightforward to verify the relations 

--^0. f--''^. T'-.,^ ,A., 

where the Jacobian of the transformation is 

J = T"f - f r') (A-5) 



d'^uj d'^uj ( d'^uj \ ' 
5^ 



2 drf' \d^drj J 



which we assume is nonzero. Note that this analysis excludes any solutions where J = 
0. The nonlinear partial differential equation (A-1) is transformed to a linear partial 



differential equation in terms of the new variables 

Q^2 Q^2 d^dr] 



^ Vt + ^S^ + 2^^7J7^ = 0- (A-6) 
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Given a solution uj{r],^) of (|A-6| ) the relations ( [A- 3D along with 

aw aw 



(A-7) 



define the solution of ( |A-1 ) paranietrically. 



We consider now solutions of the linear partial differential equation ( A-6 ) . To simplify 
the ensuing analysis we will restrict ourselves to solutions where d^Td^T does not change 
sign. This is expected to be a reasonable restriction for the vacuum structure of the theory 
( |4.1| ) since previous analysis of the actions (2.1, 2.8 , ^.15 ) suggests that d^Td'^T is either 
increasing or decreasing as T — > ztcxo. We only consider the case d^Td^T < 0, which we 
expect to most closely resemble the vacuum structure of the actions (2.1,2.8, 2.15| ). In this 
regime we can define new coordinates {p, a} by 



such that 



T] = p cosh a, i = p sinh a 



p2 = _^2 + ^2 ^ -df.Td^'T > 0, tanha = -. 



(A- 






In terms of these new variables (A-6) takes the remarkably simple form 



dp da"^ 



0. 



(A-9) 



with particular solution 



U^miV, = P""' [O^me"^" + Pme-"^"] 



for arbitrary m, Um, Pm- The most convenient way to fix the Cauchy data is to specify cv 
and dio/dp at cr = 0. In this case we write the general solution of (A-S|) as 



UJ 



{p,cr) = 7^p" (a„cosh (^^/na) + 6„sinh (^\/na)) 



(A-10) 



n=0 



where {an} determine the coefficients in the taylor expansion of w at o" = and {&«} 
determine the coefficients in the taylor expansion of duj/dp at o" = 0. 
It is straightforward to compute x{r],^), t{r],^) from ([A-3| ) and 

oo 

T{p,a) = y^p^in — 1) [an cosh (-v/ncr) + ^^sinh (-v/no")) 

n=0 

from ( [A-7| ). Because this solution is not in closed form and defined parametrically, it is 
difficult to get an intuitive feeling for the behaviour of T as a function of t and x. However, 
we believe that these general solutions do admit derivative singularities and have, for several 
choices of {a„,6„}, found points (pcCTc) corresponding to finite x and t > at which T" 
is singular though T, T' and T are regular. In this case the second derivative blows up 
because the Jacobian of the Legendre transformation ( |A-5| ) is singular. 
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